Production of 77/ Mesons in the pp — > ppr\l 
Reaction at 3.67 GeV/c 
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Abstract 

The ratio of the total exclusive production cross sections for rj/ and r\ mesons has 
been measured in the pp reaction at Pbeam = 3.67 GeV/c. The observed rj/frj ratio is 
(0.83±0.1lj^o'i|) x 10 -2 from which the exclusive r\i meson production cross section 
is determined to be (1.12±0.15^Q3 1 ) / u6. Differential cross section distributions have 
been measured. Their shape is consistent with isotropic ni meson production. 
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The study of the 77/ meson production is of particular interest because of its 
large mass compared to the other members of the ground state pseudoscalar 
meson nonet. The spontaneous breaking of chiral symmetry causes the exis- 
tence of massless Goldstone bosons, which acquire mass due to explicit chiral 
symmetry breaking, and are associated with the pseudoscalar meson nonet. 
In addition quantization effects in QCD lead to the so-called C/a(1) anomaly, 
which allows the 77/ meson to gain mass by a different mechanism than the 
Goldstone bosons [1-5]. Nevertheless, the origin of the 77/ mass and its struc- 
ture in terms of quark and gluon degrees of freedom remain controversial. 

Recent measurements of the 77/ meson by the CLEO collaboration show an 
anomalously large branching ratio of B-mesons to i]/X and rj/K [6], which 
might indicate a strong coupling of the 77/ meson to gluons [7]. Furthermore, 
the quark component of the nucleon's axial-vector matrix element measured 
in the EMC experiment [8] suggests that the rjf meson couples very weakly to 
the nucleon [9,10]. 

First measurements [11-13] of the reaction pp — > pprjf near the production 
threshold provide the possibility to determine the coupling constant g v rNN- 
However, a quantitative evaluation of this coupling constant requires answers 
to several open questions concerning the production mechanism, such as the 
roles of (i) meson-exchange currents, (ii) baryon resonances in the production 
mechanism (comparable with the role of the iV*(1535) for 7/ meson production 
[14]), and (iii) final-state interactions (FSI). 

The existing data close to threshold are consistent with different one-boson- 
exchange models (OBE) including FSI [15,16] given the ambiguities in the 
treatment of heavy-meson-exchange currents. Evaluation of the different mod- 
els requires a consistent description over a wide energy range, but data have 
been lacking at higher energies, where identification of the 77/ production can 
no longer be made solely by detecting two protons in a small forward cone [11- 
13]. 

In this letter, we report on a measurement of the 77/ / 77 production cross sec- 
tion ratio at an energy where proton-proton (pp) FSI have a much smaller 
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relative influence on the production mechanism [16,17] compared to the near 
to threshold data [11-13]. In addition, we show the differential cross sections 
of the 77/ meson as a function of the polar angle in the CM (center of mass) 
reference frame and the momentum distributions of the final state particles. 
These distributions are related to the partial wave contribution in the exit 
channel (ppr/f). 

We studied the pp reaction at the SATURNE II accelerator facility at Saclay. 
The proton beam of momentum 3.67 GeV/c was incident on a liquid hydro- 
gen target and charged products were detected using the DISTO spectrometer, 
which is described in detail elsewhere [18]. This spectrometer consisted of a 
large dipole magnet (40 cm gap size, operating at 1.0 Tm) which covered the 
target area and two sets of scintillating fiber hodoscopes. Outside the magnetic 
field, two sets of multi-wire proportional chambers (MWPC) were mounted, 
along with segmented plastic scintillator hodoscopes and water Cerenkov de- 
tectors. The scintillator hodoscopes and the Cerenkov detectors allow particle 
identification by combining either the energy loss, time of flight or Cerenkov 
light output with the particle momentum. 

The large acceptance of all detectors (~ 2° to ~ 48° horizontally and ~ ±15° 
vertically), on both sides of the beam, facilitated the coincident measurement 
of four charged particles, which was crucial for the reconstruction of many 
exclusive channels like ppK + K~ [19,20], pp7r + 7r~7r°, pKA and pKT, [21]. 

The multi-particle trigger [22], which was based on the multiplicity of the 
hodoscope elements and the scintillating fibers, selected events with at least 
three charged tracks in the final state. The results presented in this work are 
based on 4.2 x 10 7 reconstructed events with four charged particles (mainly 
ppii + n~) detected. 

The exclusive 77 meson production {pp — > ppr/) was identified via its dominant 
decay channel involving charged particles (7r + 7r~7r°, branching ratio 23.1%) 
and the reaction pp — > ppr]l was reconstructed via the decay of the 77/ meson 
into 7T + 7r~?7 (branching ratio 43.8%). The selection of the exclusive 77/ and 
77 meson production is based on two kinematical observables, the 4-particle 
(ppn + n~) missing mass (M^ iss ) and the proton-proton missing mass (M^ iss ). 
Furthermore the light output from the Cerenkov detectors together with the 
particle momentum were used to discriminate between n + mesons and protons 
in the final state. 

Since neutral pions were not detected (tt° — > 77, branching ratio 98.8%), 
pprj events were selected in which M^ iss was approximately consistent with 

a missing vr° meson (0.005 GeV 2 /c 4 < (M^ ss ) 2 < 0.035 GeV 2 /c 4 ). After 
imposing this constraint the proton-proton missing mass distribution is shown 
in Fig. 1 (upper frame). A broad signal from the 77 meson is visible near 
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M'Zss ^ 0.3 GeV 2 /c 4 . 



The assumption of a missing 7r° meson allowed a constraint ((M^ iss ) 2 = M 2 ) 
to be imposed in order to improve the (M^ iss ) 2 mass resolution by a kinemat- 
ical refit procedure. In this procedure all particle momenta were simultaneous 
re-determined under the assumption of a missing tc° meson. After applying 
the refit procedure, the resolution of the r\ meson signal in the {M^ iss ) 2 dis- 
tribution is improved by about a factor 2 (see Fig. 1, lower frame). 

In both frames the solid curve shows the sum of the signal (dotted curve) and 
the background (dashed curve). The signal line shape was taken from detailed 
Monte Carlo simulations of the detector performance and the yield was deter- 
mined by scaling the line shape to match the data using a x 2 minimization 
procedure. The small deviations at (M^ iss ) 2 ~ 0.325 GeV 2 /c 4 result from im- 
perfections of the modeling of the detector response and are included in the 
estimation of the systematic errors. 

The non-resonant reaction pp — > pp7r + 7r~7r° accounts for most of the back- 
ground under the rj meson signal. Since the exact shape of the background is 
not quantitatively known, the background has been parameterized with a 3 rd 
order polynomial, that provides the best x 2 to the fit. 

The reconstruction of the decay channel pp — ► pprjf — > ppix+ix~r\ is analogous 
to the 77 meson reconstruction described above, since the 77 meson decays 
mostly by neutral modes (branching ratio 71.5%). 

For this decay channel M^ iss and M^ iss must correspond to a missing 77/ and a 
missing 77, respectively. The (M^ iss ) 2 distribution for events with a 4-particle 
missing mass consistent with a missing 77 meson (i.e. (M^ iss ) 2 — M 2 < 

0.03 GeV 2 /c 4 ) is shown in Fig. 2 before (hatched area) and after (data points) 
the kinematical refit procedure with the constraint M^ iss = M v . The spectrum 
after the refit shows a clear signal from the 77/ meson near M 2 , ~ 0.92 GeV 2 /c 4 . 
In comparison to the reconstruction of the 77 meson the kinematical refit pro- 
cedure only improves the resolution of the 77/ signal by about 25%, due to the 
lower laboratory momenta of the outgoing protons. 

The dashed curve represents the background contribution and the solid curve 
shows the sum of the background and the signal (dotted curve). The shape 
of the signal was determined from the Monte Carlo simulations and the yield 
was determined analogously as described above for the 77 meson yield. The 
background below the 77/ meson signal originates mainly from non-resonant 
pp — > ppi]7T + 7T~ production and reactions such as pp — > pp7i + 7i^im where two 
pions are not detected. Since the exact form is unknown, the shape of the 
background was parameterized as a 4 th order polynomial, that provides the 
best x 2 to the fit. 
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The relative acceptance correction of the DISTO spectrometer for the ppr]/ vs. 
ppr] channels has been determined using Monte Carlo simulations which were 
analyzed as the measured data. 

The detector acceptance was determined as a 4-dimensional function of all 
relevant degrees of freedom in the ppr]/ and ppr] final states. 

The five degrees of freedom related to the three body decay of the 77/ and the 
77 mesons were integrated in the simulations using an isotropic orientation of 
the decay plane and the measured matrix elements [23,24]. 

After accounting for azimuthal symmetry in the production reaction, we di- 
vided the kinematically allowed phase space into four-dimensional bins and 
evaluated the efficiency for each bin separately. This was realized by storing 
the number of generated and the number of reconstructed events from the 
simulations for each phase space bin. The bin-by-bin ratio provides the effi- 
ciency correction, which was stored in the 4-dimensional acceptance correction 
matrix. 

The acceptance values were typically larger than 1 % in each phase space 
bin including branching ratios, particle decay in flight, tracking efficiency and 
particle identification efficiency The approximately 6 % efficiency loss due to 
the refit procedure was also accounted for. The data were corrected using the 
appropriate entry from this 4-dimensional acceptance correction matrix, on 
an event by event basis. For a detailed discussion of the relative acceptance 
correction method see [25]. 

The simulations indicated a very low acceptance of the apparatus for 77 mesons 
produced in the backward hemisphere in the CM (center of mass) rest frame. 
However, since the entrance consists of two identical particles a reflection 
symmetry about 0cm = 90° exists, thus the backward data are redundant for 
determining total cross sections. Therefore, we only analyzed the acceptance- 
corrected production yield in the forward hemisphere, where the acceptance 
was non-zero in each phase space bin. 

It should be noted that the acceptance correction is essentially independent of 
the event generator used in the simulations due to the complete phase space 
coverage of the DISTO spectrometer. The generator used for the correction 
assumed uniform phase space density for both reactions. This assumption was 
verified by using the same correction for different phase space populations in 
the simulations. The small deviations observed are included in the estimation 
of the systematic errors. 

Furthermore, the simulations included all important decay modes for the 77 
meson [26] for the acceptance correction of the pp — > ppr]/ — > ppix+ix~r\ and 
the pp — > ppr] reactions. Hence, background processes such as 77/ — > it + -n~r] — > 
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7r + 7r~(7r + 7r~7T°), where one or both of the observed pions are from the i] or 
events with more than four charged particles in the acceptance of the detector 
are properly accounted for. 

After correcting the 77/ meson and the rj meson production yields for the re- 
spective acceptances and branching ratios [26] as described above, the mea- 
sured total cross section ratio R = a vv ^ vvrt ,j o VP ^ vm is determined to be 
(0.83±0.1llo!il) x Where the first error is statistical and the second error 
range is due to systematic uncertainties. Because both meson channels have 
been reconstructed in events with the same four-charged-particle final state 
and measured simultaneously within the same experiment, many systematic 
uncertainties cancel when considering the production ratio. As a result, the 
systematic error is dominated by the background subtraction and the relative 
acceptance correction. 

The total cross section for the pp — > ppr\ reaction is known over a wide energy 
range [27] above and below the beam momentum of this measurement. Inter- 
polation of the existing data leads to a cross section of a^t ^ = 135 ± 35 fib 
at Pbeam = 3.67 GeV/c in good agreement of the value cr™*^, ~ 120 fib cal- 
culated by Vetter et al. [14] using a meson exchange model. By multiplying R 
by (Tpp^pp^ we obtain the total cross section o PP _ PVT] , = 1.12±0.15to;3i fib. The 
systematic error in cr^?^ is geometrically added with the systematic error 
in the pp — > pprj//pp — > pprj ratio. This result is shown in Fig. 3 (filled circle) 
together with other data closer to the production threshold [11-13] and model 
calculations (solid and dashed curves) [16]. 

The calculation from Sibirtsev et al. [16] represents a one-pion-exchange model 
including the pp FSI. The solid line denotes the full calculation and the dashed 
line corresponds to the same calculation excluding the FSI. The full calculation 
reproduces the near threshold data well, but predicts a cross section of about 
2.3 fib at our energy, which is significantly above our measurement. 

The distribution of the differential cross section as a function of the CM polar 
angle of the 77/ meson (cos(6cm)) is shown in Fig. 4. The distribution dis- 
plays no significant deviations from isotropy, indicating that the 77/ meson is 
predominantly in a s-wave state relative to the two protons. The differential 
cross section distributions were determined by producing the corresponding 
(M^ iss ) 2 spectra for each kinematical bin. Each spectrum was fitted analo- 
gously as described above to determine the yield of the signal for each bin. 
The statistical uncertainty of the yield for each bin are determined by the 
fitting procedure and are shown as vertical error bars in the differential distri- 
butions (see Fig. 4 and Fig. 5). The signal line shape was calculated for each 
spectrum individually from the Monte Carlo simulations and the background 
was also allowed to vary. 
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The relative partial wave contributions in the three particle final state pprjf 
can also be determined from the momentum distributions q and p, where q is 
the CM momentum of the 77/ meson and p is the momentum of a proton in 
the proton-proton rest frame. 

Then, the total cross section is given as the sum of the individual partial wave 
contributions [28]: 

^~E/ \M h , h \ 2 dp hh . (1) 



here the sum extends over the angular momenta h, I2 and the transition am- 
plitude for the given exit channel (M/^/J , where l\ is the orbital angular 
momentum of the two protons relative to each other and / 2 is the orbital 
angular momentum of the 77/ meson relative to the proton-proton system. 

The corresponding phase space element dpi ll2 is determined by 

d Phh ~ p 2h+1 q 2l > +2 dq, (2) 



If we assume the transition amplitude M^ 2 to be almost constant over the 
available phase space then the differential cross section as a function of the 
particle momenta are given by the variation of the phase space volume. 

The measured differential cross sections are plotted in Fig. 5. as a function of q 
(upper frame) andp (lower frame). The corresponding curves, that reproduce 
the data quite well, have been obtained from equation 2, assuming l± — l 2 — 0. 
The normalization was obtained by a simultaneous x 2 minimization to both 
differential cross section distributions. The introduction of higher values of 
I does not improve the x 2 fit to the data. The result is consistent with a 
dominant Ss-wave production of the 77/ meson, where S denotes Zi = and s 
I2 = 0, respectively. 

In conclusion, the production of the pseudoscalar 77/ meson has been studied in 
the pp reaction at pbeam = 3.67 GeV/c. The 77/ meson has been reconstructed 
by measuring its charged decay products. The 77/ / 77 cross section ratio has been 
determined and the extracted 77/ cross section has been compared to data very 
close to threshold and an one-pion exchange model including the pp FSI. While 
the model describes the data close to threshold very well it overestimates our 
data point by about 100 %. 

Furthermore, the differential cross sections indicates that the 77/ meson is pre- 
dominantly produced in a s-wave state for the two protons relative to each 
other and the 77/ meson relative to the proton-proton system. 
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Calculations up to 100 MeV above the 77/ production threshold, in the frame- 
work of a relativistic meson-exchange model [29] should be extended to higher 
energies. Comparison with our results should help to learn more about the 
different 77/ meson production mechanisms, the potential influence of a rjfp 
FS1 [30] and the magnitude of the coupling constant g vfpp . 

In this context a comparison with a recent publication on the production of 
7T°, 77 and r]l mesons in proton-proton collisions [31] should be very helpful. The 
novel approach therein factored out the pp FSI and the initial-state proton- 
proton interaction (ISI) from the total cross section to determine the phase 
space dependence of the total production amplitude for 7r° (\Aq |), rj (\Aq\) and 
rjf mesons (|Aq'|). 

Our measurement will assist to evaluate different models used for the parame- 
terization of the pp FSI which is essential for the determination of the absolute 
strength of \Aq\ 2 and hence of g v / pp - 

This paper is dedicated to the memory of Igor Falomkin. 

We acknowledge the work provided by the SATURNE II accelerator staff and 
technical support groups in delivering an excellent proton beam and assisting 
this experimental program. 
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Fig. 1. Raw spectrum of {Mj^ iss ) 2 for events with four charged particles (pp^ir^) 
in the final state and a 4-particle missing mass (M^ iss ) consistent with a missing 7r° 
meson. Before (upper frame) and after (lower frame) a kinematical refit procedure 
was applied to the data. The dashed curve represents the background contribution 
and the solid curve shows the sum of the background and of the signal (dotted 
curve) . 
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(M^ ss ) 2 [Ge W] 



Fig. 2. Raw spectrum of (M^ iss ) 2 for events with four charged particles {ppir + ir~) 
in the final state and a 4-particle missing mass {M^ iss ) consistent with a missing 
77 meson. The data points show the spectrum after a kinematical refit procedure 
was applied. The dashed curve represents the background contribution and the 
solid curve shows the sum of the background and of the signal (dotted curve). In 
addition the hatched spectrum indicates the spectrum before the refit procedure 
was applied. 
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Fig. 3. Total cross section for the pp — > ppr// reaction as a function of the available 
energy above the r]/ production threshold. Shown is the value measured in this work 
(full circle) together with other data and model calculations described in the text. 
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. 4. Differential cross section for the pp — > pprjf reaction as a function of the CM 
ar angle of the rjl meson 
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Fig. 5. Differential cross section for the pp — > ppr]/ reaction as a function of the 77/ 
meson momentum in the CM frame (q, upper frame) and the proton momentum 
in the pp-rest frame (p, lower frame). The curves denote the behavior of the three 
body phase space when the two protons are in a relative s-wave state and the 77/ 
meson is in a s-wave state relative to the protons. 
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